arXiv: 1506.01608vl [astro-ph.EP] 4Jun2015 


Cloud structure and composition of Jupiter’s troposphere from 5-pm Cassini VIMS 

spectroscopy 


R.S. Giles’"’*, L.N. Fletcher’", P.G.J. Irwin’" 

^Atmospheric, Oceanic & Planetary Physics, Department of Physics, University of Oxford, Clarendon Laboratory, Parks Road, Oxford, 

0X1 3PU, UK 


Abstract 

Jupiter’s tropospheric composition and cloud structure are studied using Cassini VIMS 4.5-5.1pm thermal emission 
spectra from the 2000-2001 flyby. We make use of both nadir and limb darkening observations on the planet’s nightside, 
and compare these with dayside observations. Although there is signihcant spatial variability in the 5-pm brightness 
temperatures, the shape of the spectra remain very similar across the planet, suggesting the presence of a spectrally-flat, 
spatially inhomogeneous cloud deck. We find that a simple cloud model consisting of a single, compact cloud is able 
to reproduce both nightside and dayside spectra, subject to the following constraints: (i) the cloud base is located at 
pressures of 1.2 bar or lower; (ii) the cloud particles are highly scattering; (iii) the cloud is sufficiently spectrally flat. 
Using this cloud model, we search for global variability in the cloud opacity and the phosphine deep volume mixing 
ratio. We find that the vast majority of the 5-pm inhomogeneity can be accounted for by variations in the thickness of 
the cloud decks, with huge differences between the cloudy zones and the relatively cloud-free belts. The relatively low 
spectral resolution of VIMS limits reliable retrievals of gaseous species, but some evidence is found for an enhancement 
in the abundance of phosphine at high latitudes. 
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1. Introduction 

The 5-pm atmospheric window is a unique region of 
Jupiter’s spectrum, where a dearth of opacity from the 
planet’s principal infrared absorbers gives us access to 
parts of the atmosphere that are otherwise hidden from 
view. The sensitivity of the 4.5-5.2 pm spectrum peaks in 
the 4-8 bar region of Jupiter’s troposphere, beneath the 
planet’s topmost cloud decks (Figured)). The observed 
brightness temperatures are therefore highly dependent on 
the properties of these clouds: the observed radiance varies 
significantly from the warm, cloud-free belts to t he cooler, 
cloudi er zones, a phenomenon first described by Westohall 
(|l969ll . This sensitivity makes the 5-pm spectral region 
extremely useful in analysing both the composition and 
cloud structure of Jupiter’s middle troposphere. 

The 5-pm data used in this study were obtained 
using t he Visual an d Infr ared Mapping Spectrometer 
(VIMS, iBrown et all 2004 1 on the Gassini spacecraft. 
In late 2000 and early 2001, Cassini performed a fly¬ 
by of Jupiter, reaching a closest approach of 9.7 million 
km (137 Jovian radii) on December 30 2000. During 
this period, VIMS made spatially-resolved spectroscopic 
measurements covering the 0.35-5.1 pm wavelength range. 
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Sromovskv and Fi^ ( 2010ll used this dataset to study fea¬ 
tures in reflected sunlight at 3 pm, but the 5-pm segment 
of this dataset has not yet been analysed. In this paper, 
we use 5-pm VIMS spectra with a range of phase and emis¬ 
sion angles to investigate Jupiter’s tropospheric composi¬ 
tion and cloud structure; we will draw conclusions about 
the cloud locations and scattering parameters, and we will 
assess the degree of constraint offered by Cassini/VIMS on 
gaseous variability. This builds on previous analyses us¬ 
ing 5-pm data from Galileo NIMS and Voyager IRIS. The 
VIMS data gives us full global coverage of Jupiter, allow¬ 
ing us to compare and contrast the 5-pm spectra from the 
belts and the zones. Additionally, VIMS made observa¬ 
tions on Jupiter’s nightside and dayside within the space 
of a few weeks, giving us the opportunity to study the 
effect of reflected sunlight. 

A typical theoretic al cloud condensation model 
from Atreva et al. ( 1997ll produces four cloud levels: NH 3 - 
ice, NH 4 SH, H 20 -ice and aqueous-ammonia solution. If a 
solar chemical composition is assumed, the bases of these 
cloud levels are located at 0.7 bar, 2.2 bar, 5.4bar and 
5.7 bar respectively. Although models such as this pro¬ 
vide valuable insights into the approximate condensation 
levels for each species, we do not necessarily expect to see 
the clouds form at these precise locations. The models 
make assumptions about both the bulk composition of the 
atmosphere and the production of photochemical prod¬ 
ucts. They also do not account for the complex dynamics 
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Figure 1; The sensitivity of the 5-(Jim window as a function 
of wavelength. The peak sensitivity is in the 4-8 bar region 
of the troposphere, beneath the planet’s topmost cloud 
decks. 


of the troposphere, neglecting the vertical mixing which 
brings up material from below and rains down material 
from above. All of these factors will act to significantly 
alter the vertical cloud structure. Observations are there¬ 
fore fundamental to understanding Jupiter’s tropospheric 
cloud decks. In the following paragraphs, we will briefly 
summarise the conclusions that have been drawn from pre¬ 
vious observational studies. 


There are many studies supporting the existence of the 
uppermost NHa-ice cloud deck. 5- and 45 -| 4 m observa¬ 
tions from Voyager IRIS were analysed bv ICierasch et al" * 


1 1986ll . who concluded that there was a highly variable 
cloud deck at ~700mbar, near the expected ammonia con¬ 
densation level. Imaging data from the Galileo space¬ 
craft showed evidenc e for thick, variable clouds in the 
750±200mbar region (IBanfield et al.L 1199811 . again consis¬ 
tent with NH 3 . Further support came in the form of spec¬ 
troscopic identification of ammonia-ice feat ures, using the 


1-3 p m region of the Galileo NIMS spectra (|BafTies_et_al 


2002) a nd 9 pm data from Gassini GIRS (|Wong et al 


2004). Baines et al. ( 2002ll found that their ammonia- 


ice spectroscopic signatures were present in less than 
1% of Jupiter’s clouds, but it has been suggested that 
this may be due to a hydrocarbon haze coating the 


ammo nia-ice particles and m asking the absorption fea¬ 
tures ( Kaloeerakis et al. . 2008[l . 

There are additional studies which support the exis¬ 
tence of both an u pper NH .3 cloud and a lower cloud 
formed of NH 4 SH. iMatcheva et al. (2005) found clouds 
in the 900-1100mbar region from Cassini GIRS obser¬ 
vations at 7.18pm, which they concluded were probably 
compose d of an upper layer o f NH . 3 and a lower layer of 
NH 4 SH. Sromovskv and F^ ( 2010ll analysed 3-pm data 
from Cassini VIMS, and found that a 500 mbar cloud com¬ 
posed of both NH 3 and NH 4 SH provided the best fit; they 
suggest that rapid upwelling could lead to the presence 
of NH 4 SH particles well above the expected condensa¬ 
tion point. Further evidence for a NH 4 SH cloud comes 


from a nalyses of the 5-pm Galileo NIMS data, lirwin et al 


( 200lh performed retrievals on the NIMS data and deter¬ 
mined that the cloud that provides the majority of the 
5 pm opaci ty is located at around 1.5 bar. This was sup¬ 
ported by Irwin and Dvudina ( 20021) . who performed a 
principal component analysis on the same data and once 
again found that the dominant opacity variations occur in 
the 1-2 bar region. Based on the location, both studies 
concluded that the dominant cloud was likely to be com¬ 
posed of ammonium hydrosulfide. 

Because of its predicted location deep in the lower tro- 
posphere, the water cloud has been difficult to observe. 


IBanfield et al.l (Il998l) used Galileo SSI observations to 
identify a deep cloud (located at a pressure greater than 
4 bar) in one region close to the Great Red Spot, which 
they c oncluded was likely t o be c omposed of water. Sim¬ 
ilarly, ISimon-Miller et ( 2000 l) found evidence for the 
cloud in very small regions of the planet by identifying a 
water ice feature near 44 pm in Voyager IRIS spectra. This 
limited detection may be due to the fact that the water 
cloud is only visible in small regions where the rest of the 
atmosphere is particularly cloud-free. 

In addition to studying the cloud structure, the 5-pm 
region allows us to retrieve abundances for tropospheric 
species and determine if there is any global variability. 
Disequlibrium species, such as PH 3 , are expected to have 
higher mixing ratios in regions of upwelling, as the gas is 
broug ht up from d e eper press ures where it is more a bun- 
dant. lirwin et~^ ( 2004 ) and Fletcher et al.l ( 20091) find 
that Gassini GIRS data shows a PH 3 enhancement at the 
equator compared to the belts on either side, supporting 
this view. Similarly, zones are expected to have a higher 
water abunda nce than belts , as m oist air is brought up 
from below; Garlson et al. ( 19921) . using Voyager IRIS 
data, found that low relative humidities were correlated 
with bright 5-pm regions, so the bright belts were found 
to be depleted in water. 

In this paper, we investigate the constraints placed on 
Jupiter’s tropospheric clouds by 5-pm VIMS data, and 
search for any spatial variability in Jupiter’s tropospheric 
gaseous species. We find that the VIMS observations can 
be modelled using a single, highly reflecting cloud layer lo¬ 
cated at a pressure of 1.2 bar or lower, and that this cloud 
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is responsible for the vast majority of the global spectral 
variability. The relatively low spectral resolution of VIMS 
limits the gaseous retrievals, but we find some evidence 
for an enhancement in the abundance phosphine at high 
latitudes. 


2. Observations 


2.1. VIMS data and calibration 

Between October 2000 and March 2001, tens of thou¬ 
sands of images of the Jovian system were taken by the 
Cassini spacecraft. These included Spm measurements 
taken by t he Visual and Infra red Mapping Spectrome¬ 
ter (VIMS, I Brown et al. . 2004). VIMS is made up of 
two imaging spectrometers, one covering the visible spec¬ 
tral range (VIMS-V) and the other covering the infrared 
(VIMS-IR). The present study uses data from the 4.5- 
5.1 pm region from VIMS-IR. 

VIMS-IR covers the wavelength range 0.85-5.1 pm with 
a spectral sampling of 16.6 nm, giving 256 wavelength 
channels. The spectral resolution varies slightly across the 
wavelength range, and has an average value of 18.7nm 
in the 5-pm window (K. Baines, personal communica¬ 
tion). The instantaneous held of view of VIMS-IR is 
0.25 X0.50 mrad. By using a two-axis scan mechanism, the 
instantaneous held of view is stepped across the scene to 
scan a full 64x64 pixel image. The total held of view is 
32 X 32 mrad, which is several times larger than the size of 
Jupiter’s disk at the distance of closest approach (Jupiter 
has a diameter of 15 mrad at distance of 9.7 million km). 

The VIMS cubes were processed using ISIS3 (Integrated 
Software for Imagers and Sp ectrometers), a software pack¬ 
age provided by the USGS ( Anderson et al.L l2004l) . T he 
VIMS Science Team calibration ( McCord et ah . 20041) is 
included in this software and geometries are assigned to 
the cubes using NASA/JPL SPICE kernels, so both radio- 
metric and geometric calibration are achieved. The hnal 
result is a data cube of dimensions 64x64x256 (one spec¬ 
tral and two spatial dimensions) and a backplane cube of 
dimensions 64x64x13, including information such as the 
latitude, longitude, phase angle and emission angle of each 
spatial point. 


2.2. Data selection 

VIMS observations were made both on the dayside of the 
planet as the spacecraft approached in late 2000 and on 
the planet’s nightside as the spacecraft moved away from 
Jupiter in early 2001. Nightside observations simply show 
the thermal emission from the planet, whereas the dayside 
observations also include a component from reflected sun¬ 
light which further complicates the analysis; Sections [Sj H] 
and [7] focus on the nightside data, and Section [5] compares 
these results to the dayside observations. 

For the nightside, a set of 30 data cubes from a 7 hour 
period on 11 January 2001 were chosen, when Cassini was 
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(a) Nightside data from 11 January 2001. 
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(b) Dayside data from 17 December 2000. 

Figure 2: Brightness temperature maps from 5.0 pm VIMS 
data on two dates, giving both nightside and dayside ob¬ 
servations of Jupiter. Smoothing has been applied to these 
maps. Black regions correspond to segments where no data 
was obtained. The CRS is located at 50°W and did not 
move significantly between the two datasets. 


a distance of 15 million km from Jupiter. The use of mul¬ 
tiple cubes vastly increases the number of available data 
points and allows us to cover more longitudes, while the 
fairly limited time period means that atmospheric condi¬ 
tions are unlikely to have changed. The date was chosen 
so as to avoid the day-night boundary (to obtain purely 
nightside emission), while maximising the spatial resolu¬ 
tion. For the subsequent dayside analysis, a set of data 
cubes from a 10 hour period on 17 December 2000 were 
used. These observations were made at a similar distance 
from the planet as the 11 January ones (16 million km), 
giving a comparable spatial resolution of ^8000 km (~ half 
the size of the CRS). 

Cylindrical projections of both the nightside and the 
dayside datasets are shown in Figure [2] In both panels, 
the difference between the warm belts and the cooler zones 
can be seen. In the nightside observations, there is a ^70 K 
brightness temperature difference between the warmest 
and the coolest regions of the planet. For the dayside ob¬ 
servations, this contrast is reduced to ^50 K; the cloudy 
zones reflect more light than the relatively cloud-free belts, 
meaning that the presence of sunlight has a greater impact 
on the cool zones than on the warm belts. Zonal mean ra¬ 
diances were then extracted from both the nightside and 
the dayside datasets. The spectra were divided into forty- 
one latitudinal bins covering the range —50° to -1-50° at 
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2.5° intervals, each with a 5° width. For each latitudinal 
bin, all spectra with an emission angle within 10° of the 
minimum value were averaged to produce one final nadir 
spectrum for each latitude segment . 

Following iFletcher et al.l ( 201lll . the final errors as¬ 
signed to these averaged VIMS spectra were 12% of the 
mean radiance in the 4.5-5.2 pm range. This is a coirserva- 
tive value, irrcluding both quadrature-estimated errors due 
to pre-flight calibration and forward-irrodel uncertaiirties 
on spectral line data. The constant error margin through¬ 
out the window prevents retrievals being weighted towards 
low-radiance regions. The size of the VIMS errors is fur¬ 
ther discussed in Section [T] 

The results of this process for the nightside data can be 
seen in Figure O Figure [5 )d shows the forty-one spectra 
from the different latitudes, all normalised to 1.0 to allow 
a comparison of spectral shapes, while Figure [3^ shows 
the absolute radiances at 5.0 pm as a function of latitude. 
While Figure shows that there is a factor of 50 differ¬ 
ence between the hottest and coldest spectra, Figure ISJd 
shows that the shape of the spectra remains almost identi¬ 
cal throughout. This cannot be explained by variations in 
the chemical composition or temperature profile of the at¬ 
mosphere. The only parameter capable of reproducing the 
huge variability in the radiance while retaining the same 
spectral shape is the cloud optical thickness, as we show 
quantitatively in the following sections. 


3. Spectral modelling 

3.1. Radiative transfer model and retrieval algorithm 

The Jupiter VIMS spectra were analysed usin g the 
NEMESIS software developed by Irwin et m] ( 200811 . This 
software has previously been used to ana l yse th e 5- 
pm VIMS spectra of Saturn ( Fletcher et al.l . I 2 OIIII and 
an e arlier version was used to analyse Jupit er NIMS 
data ( Irwin et ah . 1998 . 2001; Nixon et ah . 2001). NEME¬ 
SIS is made up of a radiative transfer code that computes 
the emergent radiation for a given atmospheric profile and 
an optimal estimation retrieval algorithm which iteratively 
determines the best-fit atmospheric parameters for an ob¬ 
served spectrum. 

The radiative transfer code solves the equation of ra¬ 
diative transfer to calculate the top-of-atmosphere radi¬ 
ation as a function of wavelength. The input informa¬ 
tion required includes a priori estimates of the atmo¬ 
spheric composition and the temperature profile, as well 
as absorption line information for each molecule. Rather 
than computing each individual absorp tion line, NEME¬ 
SIS u ses the correlated-k approximation (ILacis and Oinasl . 
1991ri to reduce the computational time. This ap¬ 
proximation allows us to reshuffle the absorption coef¬ 
ficients within a small spectral interval into ascending 
order, producing a smoothly varying funct ion that re¬ 
quire s fewer steps in numerical integration ( Irwin et al.l . 
2008[l . NEMESIS has the capability to model a multiple¬ 
scattering atmosphere. The full Mie-scattering phase 



Latitude 


(a) 5.0 Jim radiance as a function of latitude. 



(b) The normalised spectrum for each latitudinal bin. 

Eigure 3: VIMS nightside spectra, divided into 5° wide 
latitudinal bins and zonally averaged. The upper panel 
shows how the absolute radiance at 5.0 pm varies with lat¬ 
itude. The lower panel shows the full spectrum for each 
latitudinal bin, normalised to allow a comparison of their 
shapes. 
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function is approxi mated using a combin e d He nyey- was updated to have a significant ly sub-Loren t zian p rofile 


Greenstein function (|Henvev and Greensteinl . Il94lf) and 
the calc ulations are perfor med using a matrix-operator ap¬ 
proach ( Plass et 1973[l . 


was updated to nave a signincant ly sub-noren t zian p i 
in the 5-p,m window, as found bv iBaillv et al. ( 200411 . 


The NEMESIS retrieval algorithm determines the at¬ 
mospheric parameters from an observed input spectrum 
through an iterative process. Starting with an initial a 
priori atmospheric profile, synthetic “forward-models” are 
generated using the radiative transfer code. These syn¬ 
thetic spectra are compared to the observed spectrum and 
the quality of the fit is assessed by using a cost function 
comprised of two terms, the residual fit to the data and 
the deviation from the a priori atmospheric state (each 
weighted by their respective uncertainties). This cost func¬ 
tion is minimised using a Levenburg-Marquardt iterative 


scheme (Irwin et al. 


3.2. Reference atmosphere 

For this study, the Jovian atmosphere was divided into 
39 levels between 15 bar and 50mbar, equally spaced in 
log(p). The temperature and volume mixing ratio of each 
gas species is defined at each pressure level. 

The tempe rature profile i s take n from Gassini CIRS 


observations ( Fletcher et ah . 2009ll . The GIRS temper¬ 


ature profile is sensitive down to a pressure of 800 mbar, 
and has been extrapolated below this using a dry adia- 
bat, consisten t with the tempe rature profile found by the 
Galileo probe ( Seiff et id] . 1998ll (see FigureHb). The NH-^ 
and PH 3 profiles are also taken from iFletcher et ^ ( 2009[l 
and have deep volume mixing ratios of 1.862 x 10 “"^ and 
1.86 X 10“® respectively. 


The GH 4 profile was obtained from iNixon et al 


and has a deep volume mixing ratio of 1.81 x 10 The 
GH 3 D pro file assumes a constan t ratio of GH 3 D/GH 4 of 


8 X 10 ® ( Lellouch et al. . I 2 OOIII . The profile of H 2 O is 


poorly constrained. For this study, we assume a deep vol¬ 
ume mixing ratio, fixed at 1 x 10 “^ (^ the solar abun¬ 
dance) with a constant relative humidity at higher alti¬ 
tudes. For the reference profile, this relative humidity is 
set to 10%. The remaining minor atmospheric species, 
GO, GeH 4 and ASH 3 were assumed to be well mixed with 
volume mixi ng ratios of 1 . 0 ppb , 0.45 ppb and 0.24 ppb 
respectively (IBfeard et al.l . l2002n . 


3.3. Line data 

The line data sources that are used to produce the k- 
tables used in this study are described in iFletcher et al, 
( 201 lh . along with the assumed broadening parameters 
and temperature dependences. This was updated to in¬ 
clude information a bout additional GeH 4 i sotope s ob¬ 
tained using STDS ( Wenger and Ghamoionl. 1998ll. As 
with t he original GeH 4 isotope included in IFletcher et al. 
(l 201 lh . these isotopes were assumed to have a half-width 
of 0 . 1 cm~^ atm~^ and a temperat ure dependence of 


As with IFletcher et al.l (1201111 . lines were assumed to 


have a Voigt profile. The exception to this was NII 3 , which 


4. Sensitivity analysis 

The radiative transfer model described in Section o 
was used to perform a sensitivity analysis, where each pa¬ 
rameter was individually altered in order to observe its 
effect on the spectrum. Initially, a reference spectrum was 
produced using the atmospheric profile described in Sec¬ 
tion o For this sensitivity analysis, a single, spectrally- 
flat cloud layer was inserted at 0.8 bar with an optical 
thickness of 10. The atmospheric paramaters were then 
individually altered, and the new spectra were compared 
to the reference spectrum. For the molecular species, the 
entire a priori profile was scaled by factors of 0.5 and 2.0. 
For water, the deep volume mixing ratio was held constant 
at the solar abundance and the relative humidity was var¬ 
ied by factors of 0.5 and 2.0. The cloud opacity was varied 
by factors of 0.8 and 1 . 2 , and the temperature profile was 
shifted by ±20 K. The results of this process can be seen 
in Figure 01 where the initial reference spetrum is shown, 
alongside nine frames showing the difference between the 
altered profiles and the reference. 

The gaseous species with the most significant impact 
on the shape of the VIMS-resolution 5-(im spectrum are 
PH 3 and H 2 O. Although NH 3 , GH 3 D, GeH 4 , ASH 3 and 
CO all have spectral lines in this region, the relatively low 
spectral resolution means that they have a minor impact 
on the VIMS spectra. 

The single parameter that has the largest impact on the 
spectrum is the opacity of the cloud. Because the model 
cloud is spectrally flat and located well above the weight¬ 
ing function peak, the entire wavelength range is affected 
equally; changing the optical thickness scales the entire 
spectrum up or down, and a relatively minor change in 
the opacity has a significant impact on the average radi¬ 
ance. Shifting the temperature profile similarly scales the 
entire spectrum up or down, but to a lesser extent; to re¬ 
produce an effect which is comparable to a 20% change of 
in cloud opacity, the temperature shift would have to be as 
large as 20 K. A more realistic tropospheric temp erature 
differ ence 2 K between a moist and a dry adiabat (|Lewisl . 
I 995 II produces a negligible effect on the spectrum when 
compared to the variability in cloud opacity. We therefore 
fix the temperature profile to the dry adiabat used in the 
reference profile, and note that this may produce a small 
error in the retrieved cloud opacities. 

In addition to the degeneracy between cloud and tem¬ 
perature, there are further degeneracies between gaseous 
species which complicate 5-p.m retrievals at this spectral 
resolution. The regions of sensitivity for each parameter 
are not unique, but instead overlap, leading to difficulties 
in disentangling the effects of each individual parameter to 
produce a reliable retrieval. This is further illustrated by 
Figure [SI which shows the results of a retrievability test. 
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Figure 4: The sensitivity of VIMS-resolution spectra to 
changes in different atmospheric parameters. The refer¬ 
ence spectrum shows a synthetic spectrum computed for 
the atmospheric profile described in Section [XU and the 
remaining panels show the (relative) changes to this spec¬ 
trum that are caused by varying each parameter. The red 
(blue) lines correspond to a decrease (increase) of 20 K for 
the temperature, a factor of 20 % decrease (increase) for 
the cloud opacity and a factor of 2 decrease (increase) for 
the gaseous species. 


75 radiative transfer models were run with random pertur¬ 
bations applied to the reference profile. Random noise was 
then added to these synthetic spectra, in accordance with 
the 12% noise estimate for VIMS. The noisy spectra were 
then used as the input for NEMESIS retrievals. Figure [5] 
compares output was compared with the ‘true’ input val¬ 
ues, and provides the Pearson product-moment correlation 
coefficient for each parameter. A correlation coefhect of 1 
indicates a perfect positive correlation, while a coefficient 
of 0 indicates that there is no linear relationship between 
the two variables. 

The results tally with Figure 01 changes in cloud cover 
and the PH 3 volume mixing ratio have the most significant 
impact on the spectrum and can therefore be accurately 
retrieved, with Pearson product-moment correlation coef¬ 
ficients that are close to I. In contrast, the retrievals for 
the other gaseous species are poor, and the retrieved val¬ 
ues are only weakly correlated with the ‘true’ input values. 
Based on these results, the abundances of NH 3 , CH 3 D, 
GeH 4 , AsHs and CO will be held fixed at their a priori 
values during the remainder of this study. Although the 
retrievals of H 2 O are fairly poor, we have no reliable a 
priori profile for water vapour, so the relative humidity 
of water will be allowed to vary, along with the volume 
mixing ratio of PH 3 and the cloud optical thickness. 


5. Basic retrievals 


Having explored the sensitivity of the 5-pm region to 
different parameters and the potential limitations of any 
retrievals, we now turn to modelling the VIMS data itself 
using the NEMESIS retrieval algorithm. As with the sen¬ 
sitivity analysis, PH 3 and H 2 O were each allowed to vary 
via a single parameter. For PH 3 this was a scaling param¬ 
eter that effectively altered the deep volume mixing ratio, 
and for H 2 O this was the relative humidity above a fixed 
deep volume mixing ratio. For each variable, large errors 
were assigned to a priori values. Testing was carried out to 
determine whether allowing additional degrees of freedom 
in these gases would improve the fit, but this did not make 
a difference in either case. 

As an initial model, we use a single, compact, 5-km 
thick cloud, located at 0.8 bar. The cloud base pres¬ 
sure is well above the weighting function peak (4-8 bar) 
and is the approximate location of the predicted NH 3 - 
ice cloud. Figure [3] shows that the shape of the spec¬ 
tra are almost identical at all latitudes, despite large 
variations in the absolute radiance. The cloud param¬ 
eters therefore cannot vary strongly with wavelength, 
as otherwise any spectral features would be more ev¬ 
ident in the cool, optically thick spectra than in the 
warm cloud-free spectra. In this initial model, we there¬ 
fore use a spectrally flat cloud, as has been previ¬ 
ously suggested by both ground-b a sed a n d space-based 
studie s, including Drossart et al. ( I982ll. Bezard et al. 
( 1983 1. iBioraker et al.l ( 1986 1 and Roos-Serote and Irwin 
(200^. The scattering parameters for this cloud were fixed 
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Figure 5: The retrievability of different atmospheric paramaters for VIMS-resolution spectra. In each case, the true 
‘input’ value is plotted against the retrieved ‘output’ value. For H 2 O, the input and output values correspond to the 
relative humidity, for the cloud, they correspond to cloud opacity and for the remaining gaseous species, they correspond 
to the scaling factors of the a priori profile. A perfect retrieval would result in the plotted points lying along the diagonal 
line. The Pearson product-moment correlation coefficients are given in the titles. 



Figure 6 : Fit obtained using nightside data from the equa¬ 
torial zone. The VIMS data is shown in black (with the 
error shown in grey) alongside the best-ht retrieved spec¬ 
trum in red. 


at values that are broadly representative of moderately- 
sized NH 3 particles: single-scattering albedo, w = 0.9 and 
asymmetry parameter, g = 0 . 8 . 

An example of the results of these retrievals can be 
seen in Figure IH we are able to produce a good fit to 
the VIMS 5-[im data despite using a very simple cloud 
model and relatively few free parameters. The one excep¬ 
tion to the otherwise good ht is the 4.65-4.75 pm region 
of the spectrum, where there is an apparent offset be¬ 
tween the VIMS data and the retrieved spectrum. Many 
attempts were made to solve this issue, including allow¬ 
ing additional gaseous species to be retrieved, altering the 
parametrisations for the gases (including PH 3 ), inserting 
multiple cloud decks and varying the vertical cloud pro¬ 


file. Similar issues can be seen in previo us 5-pm studies of 


Jupite r usi ng Galileo/N I MS, i n cluding iRoos-Serote et al 


( 1998h and Nixon et al. ( 200lll . IPoos-Serote et all 
also had some difficulties in fitting a similar region of the 
high-resolution ISO/SWS spectra. The consistency of the 
problem suggests that the problem is with the models 
rather than the data. The mismatch is similar for both 
the cool and the warm spectra, suggesting that it is not 
due to a cloud spectral feature, and is more likely to be 
due to missing spectra l lines from a gaseous species. As 
with Roos-Serote et al.l (1998), we exclude this part of the 
spectrum to ensure that the retrievals are not driven by 
an attempt to fit this one region. 


6. Cloud constraints 

Instead of making assumptions about the composition of 
Jupiter’s clouds, we sought to determine the range of cloud 
parameters that were consistent with the VIMS data. In 
Section[5l we showed that a simple cloud model, consisting 
of a single, compact, spectrally flat cloud deck, is able to 
achieve a good fit to data. In this section, we further 
explore the cloud parameter space in order to constrain 
the cloud properties. 

6.1. Vertical structure 

Preliminary conclusions about the vertical location of 
the tropospheric cloud decks can be drawn by consider¬ 
ing the VIMS spectra alone. As described in Section [21 
Figure [3] shows that the shape of the spectra are almost 
identical at all latitudes, despite large variations in the ab¬ 
solute radiance. This immediately tells us that the main 
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Figure 7: Sensitivity of spectra to cloud locations. Each panel shows the spectrum for a variety of different cloud optical 
thicknesses. The upper panel corresponds to a compact cloud located at 0.8 bar and the lower panel corresponds to a 
compact cloud located at 2.0 bar. 



cloud decks must be located well above the 4-8 bar region 
where the weighting functions peak (see Figured]); oth¬ 
erwise, thick clouds would block out the radiation from 
some parts of spectrum more than others, giving a different 
shape to to the 5-(Jim spectrum depending on the bright¬ 
ness temperature. Additionally, having an optically thick 
cloud at a deeper, hotter part of the atmosphere will also 
start to introduce a blackbody slope into the spectrum, 
which is not observed in any of the cool spectra. These 
effects can be seen in Figure [T] which shows the results of 
inserting a cloud at two different locations in Jupiter’s tro¬ 
posphere and gradually increasing the optical thickness. In 
the upper panel, the cloud is located relatively high up in 
Jupiter’s atmosphere, and increasing the optical thickness 
decreases the observed radiance, but maintains the same 
spectral shape, just as is observed in the VIMS spectra. In 
comparison, the lower panel shows the results for a cloud 
that has been placed deeper in the atmosphere, at 2.0bar. 
Although the spectral shape remains similar for the lower 
opacities, the spectral shape starts to change for high op¬ 
tical thicknesses. In addition to changing the shape of the 
spectrum, this second case simply cannot reproduce the 
low average radiances observed in the cool zones in Fig- 
ureO at these deeper pressures, the thermal emission from 
the cloud itself is higher and even an optically thick cloud 
therefore cannot reproduce the observed low radiances. 

The constraints on the vertical location of the clouds 
were further explored by running retrievals with a wide 
range of base pressures for the single, compact cloud. This 
was done for both a cool spectrum from the EZ and a warm 
spectrum from the NEB, and the results for two of the al¬ 
titudes tested are shown in Figure [Hi When the cloud 
is located at 0.8 bar, good fits are obtained for both the 
EZ and the NEB. As the cloud is moved to progressively 
deeper pressures, the NEB hts initially remains good, but 
the EZ fits worsens significantly. The cut-off for a reason¬ 
able EZ fit is approximately 1.2 bar, while the NEB fits 
start to deteriorate for cloud base pressures of 3.0 bar or 
higher. The poor fit at deeper pressure levels is a result of 
the phenomena shown in Figured as the cloud is moved 
deeper into the atmosphere, it starts to introduce a slope 


into the spectrum that is not seen in the data. In order 
to achieve a good fit for both warm and cool spectra us¬ 
ing a single cloud deck, the clouds must be located at an 
altitude of 1.2 bar or higher. 

In addition to exploring the location of the primary 
cloud opacity, testing was also carried out to investigate 
the vertical profile of the cloud. Additional cloud decks, 
with independently variable opacities, were included in the 
retrievals, but this made very little difference to the qual¬ 
ity of the fits obtained. Similarly, extended cloud decks 
were used instead of vertically compact clouds, but this 
did not improve the fit either. 

In summary, the VIMS 5-[jim nightside data can be fit 
using a simple cloud model, consisting of a single, compact, 
spectrally flat cloud deck whose opacity varies as a func¬ 
tion of latitude, provided that this cloud deck is located 
at pressures less than 1.2 bar. It is important to note that 
this does not rule out more complicated cloud structures, 
including multiple cloud layers above 1.2 bar. The exis¬ 
tence of deeper cloud decks is also possible, provided that 
the bulk of the 5-p.m opacity still originates from the upper 
clouds. 


Similar analyses have previously been performed using 
observations from G alileo NIMS which cov e r the same 
wavelength range. Boos-Serote and IrwinI ( 2006ll used 
a grey cloud model to fit the NIMS data, and found it 
had to be located above 2.0 bar. Our study provides a 
tighter constraint (p<2.0bar) because the global VIMS 
data includes spectra from very cold regions, while the 
NIMS study was restricted to warmer regions. It is the 
cooler spectra that provide the stronger constraint on the 
cloud locat ion. Other NIMS st udies have th i s sam e dis- 
crepency; Irwin et al. (I 2 OOIII . Nixon et al. (200 ill and 


Irwin and Dvndinal (|2002ri all place their principal cloud 


decks at pressures of 1-2 bar, but since they only use 
warmer spectra, they are able to place the the clouds 
deeper in the atmosphere and still achieve good fits. 

The 5-[im VIMS results are consist ent with Cassini ob¬ 
servat ions made at other wavelengths. ISromovskv and Fry 
( 2 OIOII analysed the 3-)Jtm segment of the VIMS data. Us¬ 
ing their four-layer multiple-scattering model, they found 
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Figure 8 : Retrievals of VIMS spectra, with different assumptions about cloud location. Retrievals from the cool Equa¬ 
torial Zone are shown at the top and retrievals from the warm North Equatorial Belt are on the bottom. Each column 
corresponds to a different cloud base location. 


that the deepest, highest opacity cloud had a spatially vari¬ 
able cloud base that was located between 0.79 and 1.27 bar. 
Using a narrow spec tral window from the C IRS instrument 
centered on 7.18pm, Matcheva et al.l ( 2005 ) found that the 
cloud absorption coefficient peaks at 0.9-1.1 bar. These lo¬ 
cations are broadly consistent with the VIMS data, but we 
can neither confirm nor rule out the latitudinal variabil¬ 
ity in the cloud base pressures; we can achieve a good fit 
with a single cloud base pressure of 1.2 bar or lower, but 
a range of pressures is also possible, provided the cloud is 
located above 3.0 bar for the warm spectra from the belts 
and above 1.2 bar for the cool spectra from the zones. 


6.2. Spectral features 

As described in Section [SI the similar spectral shapes in 
Figure [3] suggest that the clouds providing the majority of 
the 5-pm opacity must be relatively spectrally flat, at least 
at the VIMS spectral resolution. Having established that 
a completely grey cloud is able to produce a good fit to the 
VIMS data, we now explore the extent to which spectral 
features are compatible which the observations. 

The two expected cloud materials in the middle tropo¬ 
sphere are NHa-ice and NH 4 SH-ice. Figure [9] shows the 
extinction cross-section and single-scattering albedo as a 
function of wavelength for a range of particle sizes. These 
functions have been calculated through Mie theory, as¬ 
suming perfectly spherical particles. From this figure, it 
is clear that the scattering and absorption parameters of 
pure particles of NH 3 and NH 4 SH vary significantly with 
wavelength within the 5-(im window. Larger particles have 
a relatively flat extinction cross-section, but a strongly 



Wavelength (|xm) 



IjimNHg - 1jimNH4SH 

lOnmNH. lOjimNH^SH 

40nni NH 3 - 40jim NH^SH 


Figure 9: Spectral parameters for NH 4 SH and NH 3 parti¬ 
cles of different sizes, calculated via Mie theory. 
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Figure 10; The effect of different cloud particles on the 
shape of the spectra. The black line in each case gives the 
fit obtained using a spectrally flat cloud. The spectral fea¬ 
tures are more pronounced in the Equatorial Zone, where 
the clouds are optically thick. 


wavelength-dependent single scattering albedo, while this 
is reversed for smaller particles. 

The effect of the wavelength-dependent cloud parame¬ 
ters on the spectral shape can be seen in Figure [TUI In 
each panel, the black line shows the model fit to the VIMS 
data that is obtained using a spectrally flat cloud. The 
remaining lines show the effect of using a “real” cloud ma¬ 
terial, rather than a grey cloud. To make the impact of the 
clouds clear, the atmospheric composition (H 2 O and PH 3 ) 
has been held constant for each cloud material, and the 
cloud opacity alone has been varied in order to produce 
an average radiance as close to the true value as possible. 
Performing a full retrieval does not however improve the 
fit sufficiently, and in some cases the algorithm attempts 
to retrieve an unphysically high abundance of the gaseous 
species in order to compensate for the poor cloud fit. 

Figure [TUI clearly shows that these “real” clouds are un¬ 
able to reproduce the shape of the VIMS spectra. This is 
particularly the case for Equatorial Zone, where the clouds 
are optically thick and any spectral features therefore 
become more prominent. The varying single-scattering 
albedo is the primary cause of these poor fits; the sharp 
gradient at ^4.65 pm seen in Figure|U]for the NH 4 SH par¬ 
ticles translates into the sharp gradient seen in EigurelTOl 
while the gradual slope for the NH 3 produces a slope in 
the spectra. 


Despite being the most plausible materials for the thick 
tropospheric clouds, pure NHs-ice and pure NH 4 SH-ice are 
not consistent with the VIMS 5-pm observations. How¬ 
ever, there are several possible explanations for this. The 
presence of several cloud decks made up of particles of 
different materials and/or sizes can act together to “blur 
out” the individual spectral features, giving a net effect of 
a roughly grey cloud. This effect, along with the fact that 
the study focussed on warmer regions where the clouds 
are thinner, allowed Irwin et "SI ( 1998ll to fit the Galileo 
NIMS 5-pm spectra using a 4-level cloud structure, made 
up of 0.5 pm tholins, 0.75 pm NH 3 particles, and 0.45 pm 
and 50 pm NH 4 SH particles. 

Alternatively, the clouds could be made of a different 
material altogether, or be masked by deposits from an¬ 
other material which causes the cloud to be spectrally flat. 
Coating of a mmonia particles by other substances has been 
suggested bv iBaines et al.l ( 2002h as an explanation for the 
absence of spectroscopically identifiable ammonia clouds in 
across the majority of the planet and iKalogerakis et aF 


1 2008h found that thin layers of hydrocarbons are able to 
alter the spectral features at 3 and 9 pm. 

The 5-pm data alone cannot distinguish between these 
various possible scenarios, but the net effect in each case 
is a roughly grey cloud. We choose the simplest option for 
reproducing this, and continue the analysis using a single 
spectrally-flat cloud. 


ff.3. Scattering properties 

The choice of scattering parameters affects the thermal 
scattering on the planet’s nightside. We now explore the 
range of scattering parameters that are consistent with the 
VIMS nightside observations. Using the nightside nadir 
data described in Section 12.21 and used in Sections 16.11 
and retrievals were performed using many values of 
(jj and g, covering the full range of physically realistic par¬ 
ticles. Although the numerical results of the retrievals var¬ 
ied with these different values (a higher single-scattering 
albedo requires a higher cloud optical thickness for the 
same spectrum), the fits produced were very uniform. 
Many different cloud parameters were capable of repro¬ 
ducing the nadir VIMS spectra, so the values of w and g 
cannot be constrained from the nadir data alone. 

Additional insights can be gained by considering obser¬ 
vations made at a range of emission angles. In the Jovian 
troposphere, temperature decreases with altitude, leading 
to limb darkening; the radiance observed decreases to¬ 
wards the edge of the disc. The extent of this limb dark¬ 
ening is heavily dependent on the cloud decks that are 
present in the troposphere; by comparing the spectra at 
different emission angles, we can constrain the scattering 
parameters of these clouds. 

The relatively low spatial resolution of the VIMS data 
means that it is not possible to isolate a single atmospheric 
feature and compare its appearance at different emission 
angles. Instead, the scatter plots in Figure [TT] show the 
limb darkening observed in the equatorial region of the 
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Figure 11: Limb darkening for different cloud parameters. 
The scatter plots show the VIMS data as a function of 
emission angle, averaged over the 5-(Jim spectral window. 
The large spread is due to the spatial inhomogeneity of the 
clouds, even within a single latitude band. Overplotted are 
the synthetic limb darkening curves for three example sets 
of cloud parameters (see Figure [12] for the full range of 
parameters tested). 
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Figure 12: The synthetic curve radiances at an angle of 
60°, for a range of scattering parameters. The red crosses 
mark the three example parameters used in Figure |TT| The 
red line marks the minimum observed radiance from the 
VIMS data; scattering parameters above this line are not 
consistent with the data. Cases where uj is low and g is 
high are ruled out. 



planet (—2.5° to 2.5°). The equatorial region was chosen 
as it provides the largest range of emission angles, and 
exhibits less spatial inhoinogeneity than other parts of the 
planet. Each point corresponds to the mean 4.5-5.2 pm 
radiance from a single pixel. 

The large spread of points is due to the inhomogeneity of 
the cloud thickness, even within a single latitudinal band 
(as seen in Figured]). Nevertheless, the general limb dark¬ 
ening trend can be seen - as the emission angles increases 
towards the limb, the radiance decreases. 

The broad scatter of points in Figure [TT] means that a 
simultaneous retrieval of data at different emission angles 
is unlikely to be productive. Two points at different posi¬ 
tions along the disc are likely to have very different cloud 
opacities, in addition to having different emission angles. 
However forward modelling does allow us to place some 
constraints on the cloud parameters that are consistent 
with the observed data. 

For a range of cloud parameter combinations (single¬ 
scattering albedo, to and asymmetry parameter, g), re¬ 
trievals were performed using nadir (low-emission angle) 
data. The retrieved parameters from the nadir data 
were then used to forward-model spectra corresponding 
to higher emission angles. The results of this process for 
three example combinations are shown by the solid lines 
in Figure im 

Each limb darkening curve is anchored to the same point 
because of the initial retrieval, but they all vary differently 
with emission angle. Despite the broad scatter of points, 
it is immediately apparent from Figure lllh that the case 
ui = 0.9, p = 0.7 is consistent with the data, but w = 
0.7, g = 0.9 is not. 

This retrieval and forward-modelling process was re¬ 
peated for more parameter combinations. For each case. 


the radiance at an emission angle of 60° was recorded and 
is plotted in Figure [12] The red line corresponds to a ra¬ 
diance of 0.062 pWcm“^ sr“^ pm“^, the lowest observed 
radiance at 60° in the VIMS equatorial data. Scattering 
parameters which produce radiances larger than this value 
are consistent with the data; scattering parameters which 
produce smaller radiances are not. Looking at Figure (TT] 
the blue curve falls within the acceptable ht region, the 
green curve is on the borderline, and the red curve does 
not give an acceptable fit. While a fairly large range of uj 
and g values give an acceptable fit to the VIMS limb dark¬ 
ening data (the bottom-right segment of the plot), this 
figure rules out the case where w is low and g is high. 

For this analysis, the location of the cloud deck was held 
constant at 0.8 bar. The effect of varying the altitude of 
the cloud deck was investigated, and it was found that no 
additional constraint was provided on top of the results 
from Section I6d1 i.e. the clouds can be located anywhere 
above 1.2 bar. The conclusions drawn from Figure IT^ are 
independent of the cloud location, provided that it is above 
the 1.2 bar level. 

Thi s is in good agre e ment with the results 
from iRoos-Serote and Irwinl ( 2006 1. who performed 
a limb darkening analysis using pairs of 5 pm data cubes 
from Galileo NIMS. With g hxed at 0.8, they found that 
an optimum fit was obtained with ui = 0.9 ± 0.5. 


7. Global retrievals 

Using the best-ht cloud model from the previous sec¬ 
tion (a single, compact grey cloud located above 1.2bar), 
retrievals were run across the entire planet in order to in¬ 
vestigate variability in the atmospheric parameters. As 
before, the free parameters in the retrievals were the cloud 
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Figure 13: North-south latitudinal retrievals of the cloud opacity, the H 2 O relative humidity and the PH 3 abundance 
from Jupiter’s nightside. The different colours correspond to the four different conditions placed on the H 2 O relative 
humidity: allowed to vary (black), fixed at 0.2% (blue), fixed at 1% (red), fixed at 3% (green). Along with the retrieved 
quantities is the goodness-of-fit as a function of latitude for each condition. 
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Figure 14: Global maps of the retrieved cloud opacity and 
PH 3 abundance from the VIMS nightside observations. 
The H 2 O relative humidity has been held fixed at 0.2%. 
Smoothing has been applied to these maps. 
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opacity, the water vapour relative humidity, and the scal¬ 
ing factor of phosphine (effectively the deep volume mixing 
ratio). Retrievals were run for each of the the forty-one 5° 
latitudinal bins described in Section and the results are 
shown in Figure IT^ The initial error estimate of 12% lead 
to underfitting the data, so the errors were reduced to 6 % 
in order to give a reduced chi-squared value of '^1 - the 
errors shown in grey are the formal errors on the retrieval, 
determined using the 6 % error value. 

7.1. Water vapour 

For the initial retrieval, the H 2 O abundance was in¬ 
cluded as a free parameter. The results for this process 
are shown in black in Figure [121 If allowed to vary, the 
zonally-averaged relative humidity ranges from ^ 0 . 2 % to 
^^3%, with a maximum around the equator and a some¬ 
what asymmetric appearance. However, Figure [5] previ¬ 
ously showed that H 2 O cannot be accurately retrieved due 
to extreme degeneracy and correlation with the other pa¬ 
rameters, so these results should be taken with caution. 

To further test these results, we then fixed the H 2 O 
relative humidity at various different levels and re-ran the 
latitudinal retrievals. Figure [T2] shows the results for three 
of these relative humidities: 0.2% (blue), 1% (red) and 3% 
(green). The fourth panel of this figure shows the impact 
that each fixed value has on the goodness-of-fit. Although 
the higher fixed relative humidities (1% and 3%) lead to 
significantly worse hts at certain latitudes, the lowest value 
( 0 . 2 %) produces hts that are comparable to the free re¬ 
trieval. We note that removing water vapour entirely from 
the retrievals signihcantly worsened the ht. 

For certain hxed H 2 O relative humidities, the cloud 
opacity and PH 3 abundance are able to adjust and pro¬ 
vide an equally good ht to the data. By considering the 
changes in the chi-square parameter, we estimate that the 
maximum hxed relative humidity allowed is approximately 
0.5%. Since we are able to achieve a good ht using a sin¬ 
gle H 2 O prohle at all latitudes (provided that the relative 
humidity is less than 0.5%) we conclude that the VIMS 
5-(im spectra do not provide any evidence for latitudinal 
variability in water vapour. Previous studies using NIMS 
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data ( R,oos-Serote et al. . 200fll) have, however, suggested 
that there is considerable local variability in the H 2 O hu¬ 
midity; since these NIMS analyses made use of high spatial 
resolution observations of small regions of the planet, this 
discrepancy may be due to sub-pixel and zonal inhomo¬ 
geneity in the relatively low resolution VIMS observations. 
Any small regions of elevated water may be rendered invis¬ 
ible by averaging over the large areas covered by the VIMS 
pixels, although a thorough exploration of gas and cloud 
degeneracies in the NIMS dataset is required to confirm 
this hypothesis. 

Assessments of spatial variability of tropospheric H 2 O 
must await higher spectral resolution space-based mea¬ 
surement to better distinguish between the competing ef¬ 
fects of water, cloud opacity and phosphine. 


1.2. Phosphine 

The latitudinal distributions of PH 3 for the different 
H 2 O conditions are shown in Figure [131 a-nd Figure [14] 
shows the global distribution for the case where the H 2 O 
relative humidity is fixed at 0.2%. These global maps 
were produced by binning the spectra into 10 ° x 10 ° size 
bins and running a retrieval at each latitude and longitude 
point. 

The different water profiles in Figure [T3] lead to the fol¬ 
lowing global averages for the PH 3 deep volume mixing 
ratio: 0.90±0.09 ppm (allowing H 2 O to vary), 0.76±0.05 
ppm (0.2% relative humidity), 0.92±0.6 ppm (1% relative 
humidity), 1.09±0.06 ppm (3% relative humidity). The 
best fits therefore suggest a deep volume mixing ratio of 
0.76-0.90 ppm. These retrieved global aver ages are consis¬ 
tent w ith previous 5-p,m studies, including iBioraker et al 


( 1986h who gives a value o f 0.7 ppm from ai rborne spec¬ 
troscopic observations, and Irwin et HI ( 1998h who gives a 
value of 0.77 ppm from Galileo NIMS. However, analyses 
of Cassini CI RS 7.7-16.6 v im ob servations have reported 
higher values: Irwin et al. ( 20041) fit the data using a deep 
volu me mixing ratio t hat va ries between I.O and 1.5 ppm, 
while Fletcher et al. ( 2009h found that values of 1.8-1.9 
ppm produced the best ht. Since CIRS is sensitive to 
higher altitudes than NIMS and VIMS, and we expect the 
phosphine abundance to decrease with altitude rather than 
increase, it is likely that this discrepancy is due to an in¬ 
consistency in the database line strengths between the two 
spectral regions. 

Although the different water vapour profiles shown 
in Figure [T3] lead to slightly different PH 3 retrievals, 
they each produce a similar latitudinal pattern: an en¬ 
hanced abundance at high latitudes compared to low lat¬ 
itudes. This is a phen omenon previously no ted in the 
northern hemisphere by Drossart et al.l ( 1990h . who used 
high-resolution S-jjim spectra to detect a 60% enhance¬ 
ment of PH 3 at high northern latitudes compared to the 
NEB; our retr i evals give an enhancement of ~60-75%. 
Drossart et al.l ( 199Clll suggest that this high-latitude en¬ 
hancement could be either due to an increase in photo¬ 



Figure 15: Zonally averaged 5.0 pm radiance as a function 
of latitude for both the nightside and the dayside. 


chemical dissociation or a decrease in the eddy diffusion 
coefficient at high latitudes. 

7.3. Clouds 

Figures [5] and [31 both showed that there was significant 
variation in 5-pm radiance with latitude; as expected. Fig¬ 
ures [TS] and [TT] shows that this is primarily due to variable 
cloud thickness. The cloud opacity retrieval shown in Fig¬ 
ure [14] is the inverse of Figure [2a] where the brightness 
temperature is high, the opacity is low, and vice versa. 
The absolute retrieved values of the cloud optical thick¬ 
ness are highly dependent on the chosen cloud parameters 
(and slightly dependent on the H 2 O relative humidity), 
but the relative changes are always reproduced, with high 
optical thicknesses in the zones (0°,±30°) and low optical 
thicknesses in the belts (±15°). 

Figure [5] also showed that there is considerable variabil¬ 
ity within latitude bands. One region that has significant 
longitudinal variation is the SEB, where the brightness 
temperature varies with distance from the Great Red Spot. 
Figure [14] shows that this is primarily due to cloud opacity. 
The thickest clouds in the SEB occur to the west of the 
GRS (50°W, 20°S), in its turbulent wake. It is not un¬ 
til the opposite side of the planet (~ 225°W), that these 
clouds thin out and the brightness temperature increases. 
This phenomenon may be explained by turbulence dredg¬ 
ing up material from deeper pressures, causing clouds to 
form. As the distance from the GRS increases, the at¬ 
mosphere becomes more quiescent and we return to the 
ordinary, relatively cloud-free appearance of a belt. 


8. Reflected sunlight analysis 

In addition to the nightside data, VIMS also made mea¬ 
surements on Jupiter’s dayside. The majority of this study 
has focussed on the nightside data, as the absence of re¬ 
flected sunlight considerably simplifies the analysis. In this 
final section, we seek to determine whether the results ob¬ 
tained from the nightside are consistent with the dayside 
observations. 

Figure [5] compared the global maps from the nightside 
and the dayside, and we noted that the contrast between 
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Figure 16: Example of the fit obtained using dayside data 
from the equatorial zone. 


the belts and the zones is smaller on the dayside than on 
the nightside. This can also be seen in Figure [TSl which 
gives the zonal averages as a function of latitude. The ad¬ 
ditional component of reflected sunlight makes little differ¬ 
ence in the warm belts and there are even points where the 
nightside is brighter than the dayside. This is because the 
belts are relatively cloud-free, so there is less reflection of 
sunlight. If we were looking at precisely the same cloud at 
different solar zenith angles, then we would expect a small 
increase due to reflected sunlight; the fact that we some¬ 
times see a small decrease is due to the fact that the gap of 
several weeks between the measurements and the low spa¬ 
tial resolution means that we are not observing identical 
atmospheric conditions. 

In the cooler, cloudier regions of the planet, such as 
the equator, the reflected sunlight component becomes 
more significant. The thicker clouds lead to more reflec¬ 
tion, and ensure that the dayside is consistently brighter 
than the nightside. At the equator, the zonal average radi¬ 
ance at 5.0pm increases from 0.25pWcm~^ sr“^ pm“^ to 


1.33pWcm“^ sr“^ pm”^, an increase of more than 500%. 
Again, part of this difference may be due to the fact that 
we are not necessarily observing identical regions of the 
planet, but reflected sunlight clearly accounts for a signif¬ 
icant part of the dayside radiance from the zones. This 
flux difference between the nightside and dayside obser¬ 
vations of_the_equatorial_zone is consistent with the anal¬ 
ysis of Drossart et al. ( 19981) . who studied the solar re¬ 
flected component of Jupiter’s 5-pm spectra from Galileo 
NIMS and found that the minimum flux level was six times 
greater on the dayside than on the nightside. 

We ran retrievals across the range of latitudes on the 
dayside, using the set of parameters from the previous 
sections: the volume mixing ratio of PH 3 , the relative 
humidity of water and the opacity of a single, compact, 
grey cloud located at 0.8 bar. For these dayside retrievals, 
the reflected solar component is included in the radiative 
transfer model. We found that we were able to obtain a 
good fit to the data using this simple cloud model; an ex¬ 
ample of the fit obtained in the equatorial region is shown 
in Figure [TBl 


The full results of the zonally-averaged retrievals are 
shown in Figure [T71 As with Figure [T31 the different 
colours refer to the different conditions placed on the H 2 O 
profile. The overall shapes of the retrievals are similar to 
the equivalent nightside plot ('Figure IT!?!): the cloud opac¬ 
ities peak at the same latitudes, and the PH 3 abundance 
still exhibits an enhancement at high latitudes. There are, 
however, a few differences. Firstly, the retrieved values of 
PH 3 exhibit additional minima which correspond to the 
peaks in cloud opacity. We believe the addition of re¬ 
flected sunlight accentuates the degeneracy between the 
cloud parameters and the PH 3 abundance, as it is in the 
PH 3 absorption wing at the short-wavelength edge of the 
window where the reflected sunlight component is most 
significant. A small change in the cloud scattering pa¬ 
rameters can lead to dramatic changes in the retrieved 
PH 3 values (a phenomenon not seen on the nightside), so 
retrievals are unreliable. Higher resolution dayside spec¬ 
troscopy would reduce this degeneracy, allowing more re¬ 
liable PH 3 retrievals. 

Secondly, although the cloud opacity has a very similar 
shape, the absolute values are slightly different. This is 
likely to be due to a combination of two factors: the opac¬ 
ities may be genuinely different, due to averaging over spa¬ 
tially inhomogeneous latitudinal band, and the retrievals 
may be mismatched due to slightly incorrect cloud scat¬ 
tering properties leading to an imbalance between thermal 
radiation and reflected sunlight. A more comprehensive 
study, taking into account shorter wavelengths, would be 
required to jointly constrain the cloud properties. 


9. Conclusions 

This paper uses the 2000/2001 observations of Jupiter 
made by the Cassini VIMS instrument in the 4.5-5.1 pm 
range to study the planet’s tropospheric composition and 
cloud structure. This builds on previous work using the 
Galileo NIMS and Voyager IRIS datasets by making use 
of (i) the full global coverage afforded by VIMS, covering 
both warm and cool regions of the planet; (ii) the combi¬ 
nation of nightside and dayside observations. We conclude 
that: 

1. VIMS 5-pm data (both nightside and dayside) can be 
modelled using a very simple cloud model, consisting 
of a single, compact, spectrally flat cloud. 

2. The bulk of the 5-pm opacity must be located suffi¬ 
ciently high in the troposphere that is does not impact 
the shape of the spectrum. For the coolest regions on 
the planet, with the thickest clouds, this requirement 
constrains the clouds to altitudes of 1.2 bar or higher. 
In warmer regions, the clouds can be placed deeper 
in the atmosphere and still achieve a good fit to the 
data. 

3. The spectra from both cloudy and relatively cloud- 
free regions have very similar spectral shapes, so the 
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Figure 17: North-south latitudinal retrievals of the cloud opacity, the H 2 O relative humidity and the PH 3 abundance 
from Jupiter’s dayside. The different colours correspond to the four different conditions placed on the H 2 O relative 
humidity: allowed to vary (black), fixed at 0.2% (blue), fixed at 1% (red), fixed at 3% (green). Along with the retrieved 
quantities is the goodness-of-fit as a function of latitude for each condition. 


clouds must be relatively spectrally-flat. Pure NH 3 - 
ice and pure NH 4 SH do not have sufficiently grey spec¬ 
tra, and are therefore inconsistent with VIMS data. It 
may be that spectral features in the clouds are masked 
by coating layers, or that multiple cloud decks act to¬ 
gether to blur out any features. 

4. The relative lack of limb darkening means that the 
cloud particles must be highly scattering. There is a 
degeneracy between the single-scattering albedo and 
the asymmetry parameter, but cases with a low single¬ 
scattering albedo and high asymmetry parameter are 
ruled out. 

5. The majority of the 5-pm global inhomogeneity can 
be accounted for by variations in the cloud opacity, 
with thick clouds in the zones and relatively cloud- 
free belts. 

6 . The retrieved globally-averaged deep volume mixing 
ratio for PH 3 was 0.76±0.05 ppm (with the H 2 O rela¬ 
tive humidity fixed at 0 . 2 %), consistent with previous 
5-(im studies. The latitudinal retrieval of PH 3 shows 
an enhancement in the abundance at high latitudes. 

7. The VIMS 5-pm spectra do not provide any evidence 
for latitudinal variability in the H 2 O relative humid¬ 
ity; if the H 2 O abundance is held fixed, the cloud 
opacity and PH 3 abundance are able to adjust in or¬ 
der to produce an equally good ht to the data. The 
low spectral resolution and high degeneracy between 
gases mean that H 2 O cannot be accurately retrieved, 
but a fixed relative humidity of less than 0.5% is able 
to provide a good fit to the data at all latitudes. 

8 . The low spectral resolution of VIMS also prevents the 
accurate retrieval of other gaseous species with weaker 
signatures in the 5-pm window: NII 3 , CH 3 D, GeH 4 , 
AsHs & CO. High resolution spectroscopy is required 
to retrieve abundances of these species and to search 
for any spatial variability. 
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